% 


o 


@ REPUBLIQUE FRAN^TMSE 

INSTITUT NATIONAL 
DE LA PROPRIETE INDUSTRIELLE 


PARIS 


(jj) N° de publication 


o 


(a n'utiliser que pour les 
commandes de reproduction) 


l) N° d'enregistrement national 


(5l) Int CI 5 : H 04 B 3/20 , H 04 N 5/60 


33 02377 


© 


DEMANDE DE BREVET D'INVENTION 


Ai 


(§) Date de depot : 10.03.93. 

r^m Prirtrito - 
r Nome . 

(Fl) Demandeur(s) : FRANCE TELECOM (Etablissement 
de droit public cree par la Loi du 2 Juillet 1990) — FR 
et TELEDiFFUSION DE FRANCE - TDF (S.A.) — FR. 

(43) Date de la mise a disposition du public de la 
w demande : 16.09.94 Bulletin 94/37. 

(7g) Inventeur(s) : Palicot Jacques, Djoko Kouam Moise et 
Veillard Jacques. 

(56) Liste des documents cites dans le rapport de 
recherche preliminaire : Se reporter a la fin du 
present fascicule. 


(60) References a d'autres documents nationaux 
apparentes : 

(73) Tituiaire(s) : 

(74) Mandataire : Ballot-Schmit. 


IN- 


UL 


(54) Procede et dispositif de filtrage d'un signal temporel numerique, et application a la correction d'echos dans 
un canal de transmission . 

(5?) Pour filtrer un signal temporel (e(t)) echantillonne par 
tmcs de N echantillons (e(n), e(k)) par une fonction de 
transfert definie dans le domaine frequentiel par LN echan- 
tillons (H(K)), le procede consiste a filtrer la fonction de 
transfert par une fenetre temporelle (g1) de largeur N puis 
a effectuer un sous-echantillonnage frequentiel de rapport 
L de facon a obtenir une fonction de transfert partielle defi- 
nie sur N echantillons (H1(k)). 

Le procede permet d'optimiser la complexite des circuits 
travaillant en temps reel. WK % 

Application notamment pour la correction des echos 
longs dans les recepteurs d'images de television. 
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PROCEDE ET DISPOSITIF DE FILTRAGE 


D 1 UN SIGNAL TEMPOREL NUMERIQUE, 
ET APPLICATION A LA CORRECTION D'ECHOS 


DANS UN CANAL DE TRANSMISSION. 


L 1 invention concerne un procede de filtrage numerique de 
signaux temporals tels que ceux re?us d'un canal de 
transmission ou d f un reseau de diffusion. 

Elle concerne egalement un dispositif de filtrage 
numerique pour la mise en oeuvre du precede. 

L' invention trouve une application particulierement 
interessante pour realiser la correction d'echos au 
niveau d'un recepteur relie a un systeme de transmission 
tel qu f un reseau de television. 

D f une fa<?on generale, un filtrage numerique concerne un 
signal temporel numerique, e'est-a-dire resultant d'un 
echantillonnage au rythme d'une certaine periode 
d 1 echantillonnage T generalement imposee en fonction des 
besoins du traitement pour tenir compte par exemple 
d'une bande passante determinee. Ainsi, par exemple, le 
signal video conforme a la norme D2-MAC est 
echantillonne avec une periode de 49,38 nanosecondes . 

Le filtrage numerique d'un tel signal peut etre effectue 
dans le domaine frequentiel, e'est-a-dire a partir de la 
trans formee de Fourier discrete du signal temporel 
numerique. Le filtrage est alors defini par une fonction 
de transfert egalement echantillonnee avec un certain 
pas d f echantillonnage en frequence df . Cette fonction de 
transfert est done la transformee de Fourier 
echantillonnee de la reponse impulsionnelle temporelle 
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du filtre a realiser. 

Mathematiquement, le filtrage consiste a calculer a 
chaque periode d f echantillonnage T le produit de 
5 convolution du signal temporel et de la reponse 
impulsionnelle du filtre. Le resultat de la convolution 
donne directement le signal filtre dans le domaine 
temporel. Si le filtre est defini par sa fonction de 
transfert dans le domaine frequentiel, la convolution se 

10 ramene au simple produit (dans 1 1 espace complexe) de la 
transformee de Fourier discrete du signal temporel a 
filtrer par la fonction de transfert du filtre pour 
chacune des valeurs de la frequence pour lesquelles la 
fonction de transfert est definie. Le resultat de ce 

15 produit complexe fournit alors la transformee de Fourier 
discrete du signal filtre, Une transformation de Fourier 
discrete inverse permet de revenir dans le domaine 
temporel. 

2 0 Pour que le filtre soit pratiquement realisable, le 
nombre d' echantillons dans les domaines temporel et 
frequentiel est evidemment limite. II en resulte des 
conditions et des limitations concernant les domaines de 
definition des grandeurs temporelles et f requentielles . 

2 5 Ainsi, conformement au theoreme de SHANNON, 

l 1 echantillonnage temporel a la periode T impose une 
limitation de la largeur de bande du signal pouvant etre 
filtre. Cette premiere condition peut s'exprimer en 
disant que la frequence d' echantillonnage 1/T d'un 

3 0 signal temporel reel doit etre super ieure au double de 

sa bande passante DF. De meme, 1 1 echantillonnage dans le 
domaine frequentiel implique une seconde condition : le 
signal temporel a filtrer n'est pris en compte que 
pendant un intervalle de temps limite DT de largeur au 
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plus egale a l 1 inverse 1/df du pas d 1 echantillonnage df 
dans le domaine frequentiel. L'intervalle maximum 1/df 
definit done une fenetre temporelle qui sera appelee par 
la suite "horizon temporel" du filtre. 

5 

La premiere des conditions pr^cedentes sera respectee si 
on limite par f iltrage la largeur de bande du signal 
temporel et/ou en choisissant une frequence 
d 1 echantillonnage 1/T suffisamment elevee. La deuxieme 
10 condition sera aussi respectee si le signal temporel 
peut etre represents par un nombre limite d 1 echantillons 
traites par le filtre a chaque periode 
d 1 echantillonnage . 

15 Si on utilise la transformation de Fourier discrete deja 
evoquee ci-dessus, un nombre limite N d 1 echantillons 
dans le domaine temporel sert a calculer un meme nombre 
N d 1 echantillons dans le domaine frequentiel* Les deux 
conditions precedentes peuvent etre alors resumees par 

20 l^s relations suivantes : 


(1) T < 1/(2 DF) 

(2) df < 1/DT 

(3) N.df > 2 DF 
25 (4) N.T > DT 

oil : 

T est la periode d' echantillonnage temporel, 
df est le pas d 1 echantillonnage frequentiel, 
N est le nombre d 1 echantillons servant pour 
3 0 la transformation de Fourier discrete, 

DF est la bande passante du signal temporel, 
DT est I 1 horizon temporel du filtre. 


En theorie, une operation de f iltrage d'un signal 


o 
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temporel au moyen d'un filtre numerique H est resumee a 
la figure 1 qui montre les operations effectuees dans 
les domaines temporel et frequentiel ainsi que les 
correspondances d'un domaine a l 1 autre. La partie 
5 super ieure de la figure fait apparaitre les divers 
parametres temporels impliques par le filtrage tandis 
que la partie inferieure montre les parametres 
correspondants (par transformation de Fourier) dans le 
domaine frequentiel. 

10 

Ainsi le signal temporel e(t) a filtrer (suppose 
continu) est echantillonne avec la periode T et convert! 
en grandeur numerique par le convertisseur A/N qui 
fournit le signal numerique echantillonne e(nT) , ou n 
15 est le rang de 1 1 echantillon correspondant e(n). Le 
filtre H etant suppose defini par sa reponse 
impulsionnelle echantillonnee h(n) , l 1 operation de 
filtrage consiste a calculer le produit de convolution : 

20 s(n) = h(n)*e(n). 

Le resultat de ce produit de convolution fournit le 
signal de sortie temporel echantillonne s(nT) du filtre. 
Une conversion numerique-analogique N/A de ce signal 
25 permet d'obtenir un signal de sortie analogique temporel 
s(t) . 

Le meme filtrage est represents dans le domaine 
frequentiel dans la partie inferieure de la figure. Le 
3 0 signal d 1 entree E(f) est le transforme du signal 
temporel e(t) . Ce signal est echantillonne avec un pas 
d 1 echantillonnage df et converti en grandeur numerique 
par le convertisseur A/N qui fournit le signal 
frequentiel echantillonne E(k.df). Le filtre H etant 
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suppose defini dans le domaine frequentiel, 1' operation 
de filtrage consiste a r4aliser le produit complexe des 
echantillons E(k.df) par les echantillons de iaeme rang k 
de la fonction de transfert echantillonnee dans le 
5 domaine frequentiel H(k.df). Les resultats de ces 
produits constituent la sortie echantillonnee dans le 
domaine frequentiel S(k.df) du filtre. Une conversion 
numerique-analogique N/A permet de passer dans le 
domaine frequentiel continu. 

10 

Comme nous I'avons deja signale precedemment, la 
realisation pratique d'un filtrage d'un signal temporel 
impose une limitation du nombre d 1 echantillons du signal 
a filtrer et de la fonction de transfert qui sont 

15 impliques dans les operations qui viennent d'etre 
decrites. La figure 2 illustre un filtrage classique 
effectue dans le domaine frequentiel oH le nombre 
d 1 echantillons est limite a la valeur N. Dans ce cas, le 
signal temporel continu e(t) se presente sous la forme 

20 d'un vecteur ou bloc d 1 echantillons [e(n)] obtenu par 
exemple par un echantillonnage a la periode T suivi 
d'une conversion analogique-numerique A/N et d'une 
conversion serie-paraliele S-P portant sur N 
echantillons successifs e(nT) • La fonction de transfert 

25 H du filtre etant elle-meme supposee echantillonnee sur 
N echantillons , elle se presente egalement sous la forme 
d f un bloc d 1 echantillons [H(k) ] . Le filtrage proprement 
dit consiste alors a effectuer en parallele les produits 
des echantillons H(k) de la fonction de transfert H par 

3 0 ceux de meme rang E(k) representant la transformee de 
Fourier discrete des signaux temporels echantillonnes 
e(n) • Cette transformee de Fourier discrete TFD opere 
sur N echantillons [e(n)] du signal temporel et peut 
etre realisee au moyen d'un circuit mettant en oeuvre un 
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algorithme du type "transformation de Fourier rapide" ou 
FFT ("Fast Fourier Transform"). Les produits S(k) des 
echantillons E(k) par H(k) definissent la transformee de 
Fourier discrete du signal temporel de sortie s(n) 
resultant du filtrage. Le signal temporel de sortie est 
obtenu par une transformation de Fourier discrete 
inverse TFD" 1 operant sur chaque bloc de N produits 
[S(k)]. Une conversion parallele-serie P-S au rythme de 
la periode d 1 echantillonnage T fournit le signal 
temporel S(t) aux instants t = nT. 

On constate que la complexity du filtre a realiser est 
directement liee a la taille des operateurs executant 
les transformations de Fourier discretes directe et 
inverse, c 1 est-a-dire au nombre N d f echantillons traites 
par ces operateurs. Par exemple, on peut demontrer que 
le nombre de multiplications necessaires pour realiser 
une transformee de Fourier rapide de taille N est egal a 
N.log2(N). Ainsi, on a done tout interet a choisir un 
nombre N le plus petit possible qui soit compatible avec 
les conditions (1) et (4) definies precedemment . II 
conviendra cependant que 1 1 echantillonnage de la 
fonction de transfert H du filtre satisfasse aussi aux 
conditions (2) et (3) . 

Dans le cas ou l'on dispose d'une fonction de transfert 
H definie dans le domaine frequentiel avec une tres 
grande precision, c 1 est-a-dire par un tres grand nombre 
d 1 echantillons H(K) separe par un pas d 1 echantillonnage 
frequentiel df tres faible, ces echantillons ne peuvent 
pas etre utilises directement car les conditions (1) a 
(4) ne sont generalement pas satisfaites. En d'autres 
termes, cela signifie que l f horizon temporel 1/df du 
filtre defini par sa fonction de transfert H(K) est 
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super ieur a I 1 horizon temporel N.T qui serait suffisant 
pour filtrer le signal. Si I 1 on suppose par exemple que 
l 1 inverse 1/df du pas d 1 echantillonnage df de la 
fonction de transfert H est un multiple L de l 1 horizon 
5 temporel NT, 1 "utilisation directe des echantillons H(K) 
necessiterait normalement de prevoir des transformations 
de Fourier discretes de taille LN, ce qui reviendrait a 
les surdimensionner inutilement. 

10 Le probleme est done le suivant. On dispose d'une 
fonction de transfert H definie sur un nombre L.N 
d 1 echantillons H(K) alors qu'il suffirait de filtrer le 
signal temporel e(n) sur N echantillons seulement. 
L 1 invention a pour but de rechercher un procede de 

15 filtrage utilisant une fonction de transfert HI definie 
sur N echantillons Hl(k) et pouvant representer avec une 
precision suffisante la fonction de transfert H, definie 
sur L.N echantillons H(K) . 

20 Dans ce but, 1' invention propose un procede de filtrage 
par une fonction de transfert d f un signal temporel 
numerique echantillonne avec une periode 

d 1 echantillonnage T et represents par sa transformee de 
Fourier discrete definie par blocs de N echantillons 

25 dans le domaine frequentiel, ladite fonction de 
transfert etant echantillonnee dans le domaine 
frequentiel et definie pour un nombre L.N 
d' echantillons, L etant un nombre entier superieur ou 
egal a 2, caracterise en ce que, prealablement au 

3 0 traitement du signal temporel, est calculee une fonction 
de transfert partielle definie dans le domaine 
frequentiel par les etapes suivantes : 

- definition d'une fonction dans le domaine 
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temporel, appelee fenetre temporelle, de valeur non 
nulle a l'interieur d'un intervalle de temps de 
duree egale a N fois la periode d 1 echantillonnage T 
et prenant une valeur nulle ou tendant vers zero a 
l'exterieur dudit intervalle, 

- calcul du produit de convolution cyclique sur L.N 
echantillons de ladite fonction de transfert 
echantillonnee par la transformee de Fourier 
discrete de ladite fenetre temporelle, 

- sous-echantillonnage dudit produit de convolution 
dans le rapport L de fagon a definir N echantillons 
de ladite fonction de transfert partielle, 

et en ce que sont effectues en temps reel les produits 
respect ivement desdits echantillons de ladite fonction 
de transfert partielle par les echantillons de meme rang 
de ladite transf ormee de Fourier discrete dudit signal 
temporel. 

On peut avantageusement definir la fenetre temporelle 
par une fonction echantillonnee dans le domaine temporel 
et dans ce cas, selon un autre aspect de I 1 invention, le 
produit de convolution cyclique consiste a effectuer un 
calcul de la transformee de Fourier discrete inverse de 
ladite fonction de transfert echantillonnee, a effectuer 
les produits respectivement des echantillons de ladite 
transformee de Fourier discrete inverse par les 
echantillons de meme rang de ladite .fenetre temporelle 
echantillonnee et a calculer la transformee de Fourier 
discrete desdits produits. 

Selon un autre aspect de 1' invention, la fenetre 
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temporelle est choisie de sorte que son integrale par 
rapport au temps soit egale a N fois la periode 
d f echantillonnage T. 

5 Le precede qui vient d'etre defini permet en fait 
d'ajuster l 1 horizon temporel du filtre a celui du signal 
temporel £ filtrer lorsque la fonction de transfert H 
initiale du filtre est trop f inement definie pour les 
besoins du filtrage. II en resulte une optimisation de 

10 la taille et de la complexity des circuits de 
transformation de Fourier discrete (tels que FFT et 
FFT" 1 ) qui doivent fonctionner en temps reel. Toujours 
dans le but de reduire la taille de ces circuits, ce 
probleme est lie a un probleme complementaire pouvant 

15 etre formule de la fagon suivante. On dispose d'une 
fonction de transfert sous la forme d' echantillons avec 
un pas d 1 echantillonnage frequentiel df donne. Pour 
filtrer un signal temporel, on souhaite utiliser des 
transformations de Fourier discretes de taille reduite N 

2 0 alors que 1' horizon temporel du signal normalement 

necessaire au filtrage est superieur a NT mais inferieur 
a 1/df. On a vu precedemment comment reduire l 1 horizon 
temporel du filtre. II reste maintenant a trouver une 
solution permettant d ! augmenter l 1 horizon temporel sans 
25 augmenter la taille N des circuits de transformation de 
Fourier discrete - 

L 1 invention a egalement pour objet un procede de 
filtrage permettant de resoudre ce probleme. Plus 

3 0 precisement, 1' invention concerne un procede de filtrage 

par ufte fonction de transfert d'un signal temporel 
numerique echantillonne avec une periode 

d 1 echantillonnage T et represents par sa transformee de 
Fourier discrete definie par blocs de N echantillons 


o 
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dans le domaine frequentiel, ladite fonction de 
transfert etant echantillonnee dans le domaine 
frequent iel et definie pour un nombre L.N 
d 1 echantillons, L etant un nombre entier super ieur ou 
egal a 2, caracterise en ce que, prealablement au 
traitement du signal temporel, sont calculees M 
fonctions de transfert partielles, definies dans le 
domaine frequentiel par les etapes suivantes : 

- definition de M intervalles de temps successifs 
ayant chacun une duree egale a N fois la periode 
d 1 echantillonnage T, 

- pour chacun desdits intervalles, definition dans 
le domaine temporel d'une fonction, appelee fenetre 
temporelle, de valeur non nulle a I 1 inter ieur de 
1 1 intervalle de temps associe et prenant une valeur 
nulle ou tendant vers zero a l f exterieur dudit 
intervalle, 

- calcul des produits de convolution cyclique sur 
L.N echantillons de ladite fonction de transfert 
echantillonnee par respectivement les transformees 
de Fourier discretes desdites fenetres temporelles, 

- sous-echantillonnage desdits produits de 
convolution dans le rapport L de fa?on a definir N 
echantillons de chacune desdites fonctions de 
transfert partielles, 

et en" ce qu'est effectuee en temps reel la somme de 
signaux frequentiels partiels organises en M blocs de N 
echantillons et resultant respectivement des produits 
des echantillons desdites M fonctions de transfert 
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partielles par des echantillons de merae rang de M blocs 
success if s de N echantillons de ladite transf ormee de 
Fourier discrete du signal temporel. 

5 Selon une variante exploitant les proprietes de 
linear ite des transf ormees de Fourier, 1' invention 
concerne egalement un procede de f iltrage par une 
fonction de transf ert d'un signal temporel numerique 
echantillonne avec une periode d f echantillonnage T et 

10 represents par sa transf ormee de Fourier discrete 
definie par blocs de N echantillons dans le domaine 
frequentiel, ladite fonction de transf ert etant 
echantillonnee dans le domaine frequentiel et definie 
pour un nombre L.N d 1 echantillons, L etant un nombre 

15 entier superieur ou egal a 2, caracterise en ce que, 
prealablement au traitement du signal temporel, sont 
calculees M fonctions de transf ert partielles, definies 
dans le domaine frequentiel par les etapes suivantes : 

2 0 - definition de M intervalles de temps successifs 

ayant chacun une duree egale a N fois la periode 
d 1 echantillonnage T, 

- pour chacun desdits intervalles, definition dans 
25 le domaine temporel d f une fonction, appelee fenetre 

temporelle, de valeur non nulle a 1 1 inter ieur de 
1 1 intervalle de temps associe et prenant une valeur 
nulle ou tendant vers zero a l'exterieur dudit 
intervalle, 


- calcul des produits de convolution cyclique sur 
L.N echantillons de ladite fonction de transf ert 
echantillonnee par respect ivement les transformees 
de Fourier discretes desdites fenetres temporelles, 


o 
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sous-echantillonnage desdits produits de 
convolution dans le rapport L de fagon a def inir N 
echantillons de chacune desdites fonctions de 
transfert partielles 

et en ce qu'est effectuee en temps reel la somme de 
signaux frequentiels partiels organises en M blocs de N 
echantillons et resultant respectivement de M produits 
successifs des echantillons desdites M fonctions de 
transfert partielles par les echantillons de meme rang 
de N echantillons de ladite transf ormee de Fourier 
discrete du signal temporel. 

L' invention permet done de decomposer une fonction de 
transfert definie sur un nombre quelconque 
d' echantillons en un banc de filtres elementaires 
definis chacun sur un nombre plus reduit d 1 echantillons . 
II en resulte que la complexity globale' des circuits 
necessaires est inferieure a celle qui resulterait de 
1' utilisation directe des echantillons de la fonction de 
transfert initiale. 

L ! invention a egalement pour objet 1 1 application du 
precede de filtrage qui vient d'etre defini pour la 
correction d'echos dans un canal de transmission. Dans 
ce cas, le filtrage porte sur un signal temporel 
numerique echantillonne present au niveau d'un recepteur 
relie audit canal de transmission, ladite fonction de 
transfert etant determinee a partir du calcul de 
1 1 inverse de la fonction de transfert dans le domaine 
frequent iel dudit canal, 

L 1 invention a egalement pour objet un filtre numerique 
pour la mise en oeuvre du procede selon 1« invention. Ce 
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filtre est caracterise en ce qu'il comporte un 
calculateur numerique programme pour calculer lesdits 
echantillons de ladite fonction de transfert partielle 
et muni d'une interface de sortie permettant de delivrer 
lesdits echantillons, et en ce que ledit filtre comporte 
des circuits cables pour calculer ladite transformee de 
Fourier discrete du signal temporel ainsi que le produit 
desdits echantillons de ladite fonction de transfert 
partielle par les echantillons de meme rang de ladite 
transf ormee de Fourier discrete du signal temporel. 

Dans le cas d'une realisation en banc de filtres, le 
filtre selon l 1 invention est caracterise en ce qu*il 
comporte un calculateur numerique programme pour 
calculer lesdits echantillons desdites fonctions de 
transfert partielles et muni d'une interface de sortie 
permettant de delivrer lesdits echantillons, et en ce 
que ledit filtre comporte des circuits cables pour 
calculer ladite transformee de Fourier discrete du 
signal temporel ainsi que les produit s desdits 
echantillons desdites fonctions de transfert partielles 
par les echantillons de meme rang de ladite transformee 
de Fourier discrete du signal temporel. 

D'autres aspects, details de realisation et avantages de 
l 1 invention apparaltront dans la suite de la description 
en reference aux figures. 

- Les figures 1 et 2 sont des schemas de principe en vue 
d'introduire les differents signaux utilises ; 


- la figure 3 est un schema de principe du procede selon 
1 1 invention ; 



2702612 


- 14 - 

- la figure 4 represente les principaux elements 
constitutifs du filtre selon 1' invention ; 

- la figure 5 represente les elements du banc de filtres 
5 selon 1 1 invention pour definir des fonctions de 

transfert partielles ; 

- la figure 6 represente des diagrammes definissant les 
fenetres temporelles utilisees dans le banc de filtres 

10 de la figure 5 ; 

- les figures 7 a 9 representent plusieurs variantes de 
realisation du banc de filtres selon 1 1 invention ; 

15 - la figure 10 represente les signaux d'horloge servant 
a synchroniser le circuit represente a la figure 9 ; 

- les figures 11 a 13 sont des schemas illustrant 
1 1 application du filtre selon l 1 invention pour la 

2 0 correction d'echos dans un canal de transmission. 

Pour faciliter la comprehension de la suite des 
explications , nous utiliserons par convention la 
notation suivante : 

25 

- L, M et N sont des entiers, avec M < L ; 

- les signaux ou fonctions de transfert (reponses 
impulsionnelles) sont en minuscules lorsqu'ils sont 

3 0 exprimes dans le domaine temporel et en majuscules dans 

le domaine frequentiel ; 

- n ou m designent le rang de 1 1 echantillon courant d'un 
signal ou d'une reponse impulsionnelle dans le domaine 
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temporel, avec o < n < N - 1 et o < m < LN - 1 ; 

- k ou K designent le rang d'un echantillon courant d'un 
signal ou d'une reponse impulsionnelle dans le domaine 
frequentiel, avec o < k < N - 1 et o < K < LN - 1 ; 

- les references entre crochets designent des blocs 
d 1 echantillons de taille N ou LN selon qu f ils renferment 
respectivement les symboles de rang n, k ou m, K. 

Les figures 1 et 2 ont servi a montrer les differentes 
grandeurs susceptibles d'etre impliquees dans un 
filtrage numerique classique* Dans ce cas, le nombre N 
d' echantillons H(k) definissant la fonction de transfert 
H du filtre est egal a celui des echantillons e(n) ou 
e(k) du signal temporel definit dans le domaine temporel 
ou frequentiel. Pour se ramener a ce cas lorsque l'on 
dispose d'un nombre super ieur LN d'echantillons H(K) de 
la fonction de transfert H definie dans le domaine 
frequentiel, on realise conformement a l 1 invention un 
traitement de ces echantillons de fa?on a definir une 
nouvelle fonction de transfert HI representee par N 
echantillons Hl(k) . Cette transformation est illustree a 
la figure 3. 

Elle consiste tout d'abord a definir dans le domaine 
temporel une fenetre temporelle gl pouvant etre 
representee par une fonction prenant la valeur 1 (a un 
coefficient d'echelle pres) a l'interieur d'un 
intervalle de temps de duree egale ^ N fois la periode 
d'echantillonnage T. On effectue alors le produit de 
convolution de la fonction de transfert echantillonnee 
[H (K) ] par la transformee de Fourier discrete [G1(K)] de 
ce tte fenetre temporelle gl de fa^on a obtenir LN 
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echantillons [H1(K)], On effectue ensuite un 
sous-echantillonnage de [H1(K)] dans le rapport L de 
fagon a obtenir la fonction de transfert desiree HI 
definie sur N echantillons Hl(k). En appliquant ces 
5 echantillons [Hl(k)] a 1 1 entree [H(k) ] du filtre 
represents a la figure 2, on obtient le filtrage desire 
du signal temporel e(t) defini sur l 1 horizon temporel 
NT. 

10 La figure 4 est un schema d* ensemble d f une variante du 
filtre mettant en oeuvre le procede precedent, Le filtre 
comporte un circuit 1 pour calculer en temps reel la 
transformee de Fourier discrete definie par blocs [E(k)] 
de N echantillons du signal temporel e(t) a filtrer. Le 

15 circuit 1 comporte un convertisseur analogique-numerique 
A/N recevant le signal temporel continu e(t) et sequence 
par un signal d'horloge CKO de periode T. Le 
convertisseur A/N fournit a chaque periode T le signal 
echantillonne e(n) applique a l 1 entree d'un 

2 0 convertisseur serie-parallele S-P synchronise par un 

signal d'horloge CK de periode NT de fagon a fournir des 
blocs de N echantillons [e(n)]« En traitant ce bloc par 
un circuit de transformation de Fourier discrete FFT-N 
de dimension N, on obtient la transformee de Fourier 
25 discrete [E(k)]. Bien entendu, le circuit 1 serait 
superflu si l'on disposait par d'autres moyens des 
echantillons E(k) . 

La figure 4 montre egalement une fagon particuliere de 

3 0 realiser la convolution de la fonction de transfert 

initiale [H (K) ] par la fenetre temporelle gl definie 
precedemment . Selon cette realisation, la fenetre 
temporelle gl est sous la forme d'une fonction 
echantillonnee gl(m) definie sur LN echantillons dans le 
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domaine temporel. On calcule alors par une transformee 
de Fourier inverse la reponse impulsionnelle [h(m)] de 
[H(K) ] . On peut utiliser par exemple un operateur du 
-type "transformee de Fourier rapide inverse" FFT -1 -LN de 
dimension LN. Cette transformation fournit LN 
echantillons [h(m)] qui sont multiplies par les 
echantillons de meme rang de la fenetre temporelle 
echantillonnee [gl(m)]. Une transformation de Fourier 
FFT-LN de ces produits Jhl(m)] fournit LN echantillons 
[H1(K)]. Comme precedemment, ces Echantillons sont 
sous-echantillonnes dans un rapport L de fa<?on a fournir 
N echantillons [Hl(k)]. 

Le circuit comporte des multiplieurs pour effectuer en 
parallele et en temps reel les produits des echantillons 
Hl(k) par les echantillons de meme rang k de la 
transformee de Fourier discrete E(k) du signal temporel. 
Ces produits fournissent les echantillons [Sl(k)] qui 
constituent la transformee de Fourier discrete du signal 
de sortie du filtre. A partir des echantillons [Sl(k)], 
on peut obtenir le signal de sortie echantillonne dans 
le domaine temporel s(n) au moyen d'une transformation 
de Fourier discrete inverse FFT""i-N. 

La figure 5 represente une generalisation du filtre 
precedent sous la forme d'un banc de filtres realisant 
une decomposition de la fonction de transfert initiale 
H(K) en M fonctions de transfert partielles HI, H4, 
... definies chacune sur N echantillons, avec M < L. On 
definit pour cela M fenetres temporelles gl, g2, g^ 
telles que representees sur le diagramme de la figure 6. 
Chacune de ces fenetres correspond a un horizon temporel 
de largeur NT et la juxtaposition de M fenetres 
adjacentes permet done de definir un horizon temporel de 
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largeur M.NT. II est a noter que les fenetres 
temporelles ne sont pas necessairement des fenetres 
rectangulaires mais peuvent avoir 1' allure representee 
en pointilles sur le diagramme de fagon a deborder de 
5 part et d 1 autre de la fenetre rectangulaire de reference 
tout en conservant la meme surface. Cette disposition 
permet en particulier d'eviter des difficultes lors de 
la recomposition du signal de sortie a partir des 
signaux de sortie partiels obtenus par le filtrage du 
10 signal temporel d 1 entree par les fonctions de transfert 
partielles HI, H2, . .., H^- 

Pour chacun des etages du banc de f litres, on effectue 
les operations qui ont ete decrites en reference aux 

15 figures 3 et 4 . Ainsi, selon une premiere possibility, 
on calcule pour chaque etage le produit de convolution 
cyclique de la fonction de transfert H(K) par la 
transformee de Fourier discrete G1(K) , . . . , G4 (K) , 
de la fenetre temporelle gl, . .., g4, ... de l 1 etage 

2 0 considere. Selon une autre possibility, le produit de 
convolution est obtenu par le calcul de la transformee 
de Fourier discrete inverse h(m) de la fonction de 
transfert H(K) suivi du calcul des produits hl(m), 
h4 (m) , ... respectivement des echantillons h(m) par les 

2 5 echantillons de meme rang de chaque fenetre temporelle 
echantillonnee gl(m), g4(m) / ... definie dans le 
domaine temporel. Un calcul de la transformee de Fourier 
discrete de ces produits hl(m) , . . . , h4 (m) , ... fournit 
les blocs d 1 echantillons [H1(K)], , [H4 (K) ] , ... . 

30 

Comme" precedemment , les echantillons H1(K) , . . . , H4 (K) , 
. . . subissent un sous-echantillonnage de rapport L pour 
fournir les fonctions de transfert partielles Hl(k) , 
. . . , H4 (k) , ... recherchees. 
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Les fonctions de transfert partielles HI, H4, . 

peuvent etre utilisees pour reconstituer un signal de 
sortie defini par blocs de N echantillons seulement mais 
prenant en compte un horizon temporel super ieur a NT, La 
fa?on de realiser cette recomposition par une cellule de 
f iltrage va maintenant etre decrite en reference aux 
figures 7 a 10. 

Selon une premiere possibility representee a la figure 
7, la cellule de f iltrage 2 constitue M blocs successifs 
de N echantillons de la transformee de Fourier discrete 
[E(k)] au moyen de circuits a retard montes par exemple 
en registre decalage. Ces circuits a retard fournissent 
en sortie des blocs decales les uns par rapport aux 
autres d'un intervalle de temps egal a NT. M 
multiplieurs effectuent les produits des echantillons 
des fonctions de transfert partielles [Hl(lc) ] , • . . , 
[H4(k)], ... par les echantillons de meme rang des blocs 
successifs issus des circuits a retard. Les multiplieurs 
fournissent en sortie M blocs de signaux frequentiels 
partiels de N echantillons [Sl(k)] f [S4(k)] 7 ... qui 
sont ensuite sommes en parallele dans un additionneur 
dont la sortie fournit N echantillons de la transformee 
de Fourier discrete S(k) du signal de sortie du filtre. 

II est important de noter que si l'on prend M = L, on 
realise un f iltrage sur un horizon temporel du signal 
d 1 entree egal a l 1 horizon temporel de la fonction de 
transfert initiale H(K) . 

La figure 8 represente une seconde variante de 
realisation de la cellule de f iltrage 2, 
f onctionnellement equivalente a celle de la figure 7. 
Ces deux realisations se distinguent I 1 une de 1' autre 


o 


2702612 


- 20 - 

par le fait que les produits et les retards sont 
permutes . 

Selon une autre variante representee a la figure 9 et 
correspondant a titre d'exemple a un banc de filtres a 
quatre etages, les operations de retard et de soiwnation 
sont alternees de fagon a remplacer la sommation de M 
blocs d 1 echantillons par une succession de sommations de 
deux blocs seulement. Ceci presente 1 1 avantage de 
simplifier la realisation des somateurs et par 
consequent de reduire le temps de calcul. 

Le f onctionnement du circuit de la figure 9 se deduit 
immediatement du schema et du chronogramme de la figure 
10 qui montre le signal d'horloge d 1 echantillonnage CKO 
et le signal d'horloge de blocs CK utilise pour 
synchroniser le circuit. 

Pour illustrer une application particulierement 
interessante du procede de filtrage selon 1' invention, 
la figure 11 represente un schema general d'un exemple 
de systeme de transmission de donnees, tel qu'un systeme 
de transmission d 1 images de television. 

De fagon classique, le systeme est constitue d ! un 
emetteur 3, d'un reseau de transmission 4 et d f un 
recepteur 5. Le recepteur 5 est constitue du montage en 
cascade d'un demodulateur 6, d'un filtre passe-bas 7, 
d'un convertisseur analogique-numerique 8 associe a un 
egaliseur 9 et eventuellement d'un filtre H de 
correction d'echos avantageusement realise conformement 
a I 1 invention. L' egaliseur 9 est prevu a la fois pour la 
recuperation du signal d'horloge du signal regu et pour 
realiser une suppression des echos dits "courts" qui 
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sont inevitablement produits dans 1" ensemble du canal de 
transmission Q. La suppression des echos courts par 
l'egaliseur 9 a cependant pour contrepartie la creation 
d 1 echos dits "longs" qui bien que tres attenues nuisent 
a la qualite de 1" image. C'est la raison pour laquelle 
il est utile de prevoir le f iltre supplementaire H 
destine a supprimer ces echos longs. Cependant le 
traitement des echos longs suppose par definition que le 
f iltre H devra prendre en compte un horizon temporel 
important du signal a filtrer. C'est pourquoi le procede 
et le f iltre selon l 1 invention sont particulierement 
bien adaptes pour resoudre ce probleme sans necessite de 
prevoir des filtres volumineux. 

La figure 12 represente une realisation d' ensemble du 
f iltre H utilisable notamment pour la correction d 1 echos 
longs. 11 est essentiellement compose d'une cellule de 
filtrage 2 conforme a l 1 invention relie d'une part a un 
calculateur 10 et au signal d 1 entree [e(n) ] par 
1 ' intermediaire d'un circuit de transformation de 
Fourier rapide FFT-N. Le calculateur 10 regoit le signal 
d 1 entree en vue de calculer la fonction de transfert Q 
du canal a partir d'un signal de wobulation transmis par 
l'emetteur 3. Le signal de wobulation permet au 
calculateur de calculer avec une bonne precision la 
fonction de transfert echantillonnee dans le domaine 
frequent iel du canal Q. Le calculateur peut alors 
calculer l 1 inverse de cette fonction Q pour obtenir la 
fonction de transfert du filtre H apres une correction 
eventuelle en vue d 1 assurer la stabilite du filtre. 
Avant'ageusement, le calculateur 10 pourra mettre en 
oeuvre un algorithme d 1 adaptation permettant d'ajuster 
regulierement les coefficients de la fonction de 
transfert calculee. On utilisera par exemple la ligne 
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624 ciu signal D2-MAC pour transmettre le signal de 
wobulation. 


Enfin, le calculateur sera programme pour calculer les 
5 coefficients des fonctions de transfert partielles HI, 
• . . , H4, . . . conformement au procede selon l 1 invention. 
Bien entendu, on pourra utiliser le meme calculateur 
avec une cellule 2 reduite a un multiplieur si une seule 
fonction de transfert HI suffit. Dans le cas contraire, 
10 ii faudra prevoir M fonctions de transfert partielles et 
une cellule 2 prevue en consequence de sorte que le 
produit M.N, T. soit au moins egal a la longueur de 
1 1 echo a corr iger . 

15 A titre d'exemple de realisation, la figure 13 montre 
les principaux elements constitutifs du calculateur 10. 
Le calculateur est organise autour d'un bus B sur lequel 
sont connectes un processeur 11, une memoir e de 
programme 12, une premiere memoire 13 pour contenir les 

2 0 echantillons WOB(n) du signal de wobulation et une 

seconde memoire 14 prevue pour contenir les coefficients 
calcules du filtre H ainsi que ceux des fonctions de 
transfert partielles Hi. Une interface d 1 entree IE et 
une interface de sortie IS permettent au bus B de 
25 communiquer avec respectivement le signal d 1 entree 
[e(n)j et les multiplieurs de la cellule de filtrage 2. 

Le calculateur qui vient d'etre deer it etant de type 
classique, sa constitution detaillee et son 

3 0 f onctionnement ne necessitent pas d 1 explications 

supplementaires compte tenu des indications deja donnees 
pr ecedemment . 

Cette realisation mixte qui combine un calculateur 
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programme 10 pour le calcul des coefficients et des 
circuits cables pour le traitement en temps reel du 
signal vise a optimiser les coats de realisation tout en 
satisfaisant aux criteres de performance. 
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30 
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R E VEND I CAT I ONS 

1. Precede de filtrage par une fonction de transfert (H) 
d'un signal temporel numerique ([e(n)]) echantillonne 
5 avec une periode d" echantillonnage T et represents par 
sa transformee de Fourier discrete ([E(k)]) definie par 
blocs de N echantillons dans le domaine frequentiel, 
ladite fonction de transfert (H) etant echantillonnee 
( [H (K) ] ) dans le domaine frequentiel et definie pour un 

10 nombre L.N d 1 echantillons, L etant un nombre entier 
superieur ou egal a 2, caracterise en ce que, 
prealablement au traitement du signal temporel ([e(n)]), 
est calculee une fonction de transfert partielle 
([Hl(k)]) definie dans le domaine frequentiel par les 

15 etapes suivantes : 

definition d'une fonction dans le domaine 
temporel, appelee fenetre temporelle (gl) , de 
valeur non nulle a l 1 inter ieur d'un intervalle de 

2 0 temps de duree egale a N fois la periode 

d 1 echantillonnage T et prenant une valeur nulle ou 
tendant vers zero a l'exterieur dudit intervalle, 

calcul du produit de convolution cyclique 
25 ([H1(K)]) sur L.N echantillons de ladite fonction 

de transfert echantillonnee ( [H(K) ] ) par la 
transformee de Fourier discrete ([G1(K)]) de ladite 
fenetre temporelle (gl) , 

3 0 - sous-echantillonnage dudit produit de convolution 

X[H1(K)]) dans le rapport L de fa?on a definir N 
echantillons de ladite fonction de transfert 
partielle ( [HI (k) ] ) , 
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et en ce que sont ef fectues en temps reel les produits 
([Sl(k)]) respectivement desdits echantillons de ladite 
fonction de transfert partielle ([Hl(k)]) par les 
echantillons de meme rang de ladite transf ormee de 
Fourier discrete ([E(k)]) dudit signal temporel 
(Ce(n)]) . 

2. Procede de filtrage selon la revendication 1, 
caracterise en ce que ladite fenetre temporelle (gl) 
etant une fonction echantillonnee, ledit produit de 
convolution cyclique ([Hl(k)]) consiste a effectuer un 
calcul de la transf ormee de Fourier discrete inverse 

([h(m)]) de ladite fonction de transfert echantillonnee 

( [H(K) ] ) , a effectuer les produits ( [hi (m) ] ) 
respectivement des echantillons de ladite transf ormee de 
Fourier discrete inverse ([h(m)]) par les echantillons 
de meme rang de ladite fenetre temporelle echantillonnee 

(gl) et a calculer la transf ormee de Fourier discrete 

([H1(K)]) desdits produits ([hl(m)]). 

3. Procede de filtrage selon l f une des revendications 1 
ou 2 caracterise en ce que ladite fenetre temporelle 
(gl) est choisie de sorte que son integrale par rapport 
au temps soit egale a N fois la periode 
d'echantillonnage T. 

4. Procede de filtrage par une fonction de transfert (H) 
d'un signal temporel numerique ([e(n)]) echantillonne 
avec une periode d'echantillonnage T et represents par 
sa transf ormee de Fourier discrete ([E(k)]) definie par 
blocs' de N echantillons dans le domaine frequentiel f 
ladite fonction de transfert (H) etant echantillonnee 
( [H(K) ] ) dans le domaine frequentiel et definie pour un 
nombre L.N d 1 echantillons, L etant un nombre entier 
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super ieur ou egal a 2, caracterise en ce que, 
prealablement au traitement du signal temporel ([e(n)]>, 
sont calculees M fonctions de transfert partielles 

([Hl(k)], , [H4(k)]) definies dans le domaine 

frequentiel par les etapes suivantes : 

- definition de M intervalles de temps successifs 
ayant chacun une duree egale a N fois la periode 
d 1 echantillonnage T , 

- pour chacun desdits intervalles, definition dans 
le domaine temporel d'une fonction, appelee fenetre 
temporelle (gl f . . . , g4) , de valeur non nulle a 
l'interieur de 1 1 intervalle de temps associe et 
prenant une valeur nulle ou tendant vers zero a 
l'exterieur dudit intervalle, 

calcul des produits de convolution cyclique 
([H1(K), H4(K)]) sur L.N echantillons de 

ladite fonction de transfert echantillonnee 
([H(K)]) par respectivement les transformees de 
Fourier discretes ( [Gl (K) , . • . , G4 (K) ] ) desdites 
fenetres temporelles (gl, g4) , 

sous-echantillonnage desdits produits de 
convolution ([H1(K), H4 (K) ] ) dans le rapport L 

de fagon a definir (gl, g4) N echantillons de 

chacune desdites fonctions de transfert partielles 
([Hl(k) , • . . , H4(k) ]) 

et en'" ce qu'est effectuee en temps reel la somme de 
signaux frequentiels partiels organises en M blocs de N 
echantillons ([Sl(k), . S4(k)]) et resultant 
respectivement des produits des echantillons desdites M 
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fonctions de transfert partielles ([Hl(k), . , H4(k)]) 
par des echantillons de meme rang de M blocs successifs 
de N echantillons de ladite transf ormee de Fourier 
discrete ([E(k)]) du signal temporel ([e(n)]). 

5. Procede de filtrage par une fonction de transfert (H) 
d f un signal temporel numerique ([e(n)]) echantillonne 
avec une periode d 1 echantillonnage T et represents par 
sa transf ormee de Fourier discrete ([E(k)]) definie par 
blocs de N echantillons dans le domaine frequentiel, 
ladite fonction de transfert (H) etant echantillonnee 
( [H(K) ] ) dans le domaine frequentiel et definie pour un 
nombre L.N d* echantillons , L etant un nombre -entier 
super ieur ou egal a 2, caracterise en ce que, 
prealablement au traitement du signal temporel ([e(n)]), 
sont calculees M fonctions de transfert partielles 
( [HI (k) ] , ... r [H4 (k) ] ) def inies dans le domaine 
frequentiel par les etapes suivantes : 

- definition de M intervalles de temps successifs 
ayant chacun une duree egale a N fois la periode 
d 1 echantillonnage T, 

- pour chacun desdits intervalles, definition dans 
le domaine temporel d f une fonction, appelee fenetre 
temporelle (gl, g4) , de valeur non nulle a 
I'interieur de l f intervalle de temps associe et 
prenant une valeur nulle ou tendant vers zero a 
l f exterieur dudit intervalle, 

- calcul des produits de convolution cyclique 
( [HI (K) , • , . , H4 (K) ] ) sur L.N echantillons de 
ladite fonction de transfert echantillonnee 
( [H(K) ] ) par respectivement les transf ormees de 
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Fourier discretes ([G1(K), , G4 (K) ] ) desdites 

fenetres temporelles (gl, . g4) , 

sous-echantillonnage desdits produits de 

convolution ([H1(K), , H4 (K) ] ) dans le rapport L 

de fagon a definir (gl, , g4) N echantillons de 

chacune desdites fonctions de transfert partielles 
([Hl(k), . .., H4(k)]) 

et en ce qu'est effectuee en temps reel la somme de 
signaux frequentiels partiels organises en M blocs de N 
echantillons ( [Sl(k) , . • , , S4 (k) ] ) et resultant 
respectivement de M produits successifs des echantillons 
desdites M fonctions de transfert partielles 
([Hl(k), H4(k)]) par les echantillons de meme rang 

de N echantillons de ladite transf ormee de Fourier 
discrete ([E(k)]) du signal temporel ([e(n)]). 

6. Procede de filtrage selon l'une des revendications 4 
ou 5, caracterise en ce que lesdites fenetres 
temporelles (gl, g4) etant des fonctions 

echantillonnees r chacun desdits produits de convolution 
cyclique ([Hl(k)]) consiste a effectuer un calcul de la 
transf ormee de Fourier discrete inverse ([h(m)]) de 
ladite fonction de transfert echantillonnee ([H(K)]) / a 
effectuer les produits ([hl(m)]) respectivement des 
echantillons de ladite transformee de Fourier discrete 
inverse ([h(m)]) par les echantillons de meme rang de 
chacune desdites fenetres temporelles echantillonnees 

(gl) et a calculer la transformee de Fourier discrete 

([H1(K)]) desdits produits ([hl(m)])- 

7o Procede de filtrage selon l'une des revendications 4 
a 6 , caracterise en ce que chacune desdites fenetres 
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temporelles est choisie de fagon a ce que son integrale 
par rapport au temps soit egale a N fois la periode 
d'echantillonnage T et de sorte que la somme de deux 
fenetres adjacentes soit egale a l 1 unite, 

8. Application du procede de filtrage selon l'une des 
revendications 1 a 7 pour la correction d'echos dans un 
canal de transmission, caracterisee en ce que le 
filtrage porte sur un signal temporel numerique 
echantillonne present au niveau d'un recepteur relie 
audit canal de transmission, ladite fonction de 
transfert (H) etant determinee a partir du calcul de 
l f inverse de la fonction de transfert (Q) dans le 
domaine frequentiel dudit canal. 

9. Application du procede de filtrage selon la 
revendication 8, caracterisee en ce que ladite fonction 
de transfert du canal est determinee au moyen d'un 
signal de wobulation produit par un emetteur relie a 
l 1 autre extremite dudit canal. 

10. Application du procede de filtrage selon l'une des 
revendications 4 a 7 pour la correction d'echos dans un 
canal de transmission, caracterisee en ce que ledit 
signal temporel numerique ([e(n)]) est un signal present 
au niveau d'un recepteur relie a une extremite dudit 
canal, en ce que ladite fonction de transfert (H) est 
1' inverse de la fonction de transfert (Q) dudit canal, 
en ce que lesdits M intervalles de temps sont adjacents 
et en ce que les nombres M et N sont choisis de f agon a 
ce que le produit M.N. T soit au mo ins egal a la longueur 
de I 1 echo a corriger. 

11. Filtre numerique pour la mise en oeuvre du procede 


o 


o 


2702612 


- 30 - 


de filtrage selon I'une des revendications 1 a 3, 
caracterise en ce qu'il comporte un calculateur 
numerique programme pour calculer lesdits echantillons 
de ladite fonction de transfert partielle ([Hl(k)]) et 
muni d f une interface de sortie permettant de delivrer 
lesdits echantillons, et en ce que ledit filtre comporte 
des circuits cables pour calculer ladite transformee de 
Fourier discrete ([E(k)]) du signal temporel ([e(n)]) 
ainsi que le produit desdits echantillons de ladite 
fonction de transfert partielle ([Hl(k)]) par les 
echantillons de meme rang de ladite transformee de 
Fourier discrete ([E(k)]) du signal temporel ([e(n)]). 

12. Filtre numerique pour la mise en oeuvre du procede 
de filtrage selon l'une des revendications 4 a 7 , 
caracterise en ce qu'il comporte un calculateur 
numerique programme pour calculer lesdits echantillons 
desdites fonctions de transfert partielles ([Hl(k)]) et 
muni d'une interface de sortie permettant de delivrer 
lesdits echantillons, et en ce que ledit filtre comporte 
des circuits cables pour calculer ladite transformee de 
Fourier discrete ([E(k)]) du signal temporel ([e(n)]) 
ainsi que les produits desdits echantillons desdites 
fonctions de transfert partielles ([Hl(k)]) par les 
echantillons de meme rang de ladite transformee de 
Fourier discrete ([E(k)]) du signal temporel ([e(n)]). 
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